Casimir Effect in Fluids above the Isotropic-Lamellar Transition 
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We study fluctuation-induced interaction in confined fluids above the isotropic-Iamellar transition. 
At an ideally continuous transition, the disjoining pressure has the asymptotic form n(d oo) ~ 
—CkBTq^/d, where d is the interwall distance, qo is the wavenumber of the scattering peak, and 
C — l/(47r) in the strong anchoring limit. The long-rangedness is enhanced due to continuous 
distribution of soft modes in the q-space. An unconventionally strong Casimir force with a range of 
several lamella thicknesses is realistic above the transition. We also flnd an oscillatory force profile 
near a surface-induced transition. 

PACS numbers; 68.15.-He, 64.60.Cn, 61.25.Hq, 82.70.-y 
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Fluids in confined geometries show a variety of 
phenomena that are not observed in bulk. Among 
them, thermal-fluctuation-induced interactions between 
boundary walls or objects immersed in the fluid have 
attracted much attention. They are called the Casimir 
forces by analogy with the quantum orig inal 0. The 
forces are long-ranged in systems with soft modes, such 
as critical simple fluids [ 2|,|j and liquid crystals , with 
the decay law and amplitude modulate d by various sur- 
face effects . Recent evidence |^,§,|0[ supports the 
view that the interaction is ubiquitous in correlated flu- 
ids 0. 

In this Letter, we address the effect in structured fluids 
above the isotropic-Iamellar (I-L) transition. The model 
system we consider has two characteristic lengths, which 
describe decay and oscillation of the correlation func- 
tion. Physical realizations of the model include block 
copolymers in the disordered phase and bicontinuous mi- 
croemulsions. Films of block copolymers have consti- 
tuted the subject of numerous papers [ p^p^|Jl^JT5| , p^ . 
Most of them focus on the lamellar phase, while a few 
theoretical works treat the system above the I-L tran- 
sition The latter works discuss the effect of a 
surface field, which induces a mean-field interaction be- 
tween boundary walls. In comparison to the case of 
copolymers, much less is known about confined fluids 
containing short-chain surfactants [p"7|jl^jl9|] . A mean- 
field study with a standard Ginzburg-Landau model of 
microemulsions found a surface-induced I-L transition 
that preempts the bulk transition p7[ . 

The bulk property of the system we assume is described 
by the Hamiltonian 



(1) 



where e is a constant, is the order parameter, and 
Po and Qq are characteristic wavenumbers that gener- 
ally depend on the temperature. This and similar types 
of Hamiltonians, with or without additional nonlinear 
terms, have been applied to various kinds of complex 
fluids [|0|j2|j2|,||]. For symmetric AB-block copoly- 
mer melts, for instance, represents the excess of A- 



monomer's (say) volume fraction with respect to its 
spatial average, and the parameters are given by e = 
0m27kBTRg/N^/^, qo = 1.95/Rg, and po = 1.43[(xc - 
XjNY/'^ / Rg. Here, Rg is the chain's gyration radius, N 
is the polymerization index, x is the Flory interaction 
parameter, and Xc = 10.5/Af locates the mean-fleld TL 
transition ||2C| ]. To be general, we shall use the ratio po/9o 
as the measure of distance from the transition tempera- 
ture. 

We consider a confined film of thickness d and of 
macroscopic area. The total Hamiltonian is the sum of 
bulk and surface contributions, Ti = Hb + Ti-s- We write 
the bulk Hamiltonian (Q) in the real space as 



nb = - dr 

^ JO<z<d 



(VW~2g2(v</>)2 + (p4 + 54)^2 



(2) 



In this form, the Hamiltonian is equivalent to the 
Teubner-Strey model of microemulsions [pl| in the tem- 
perature region between the bicontinuous-lamellar tran- 
sition and the Lifshitz point. For chemically neutral and 
identical walls, the surface part of the Hamiltonian can 
be expanded up to the second order in (j) and V0 as EJ] 



ujs(p^ + 5.(V0)2] [S{z) + S{z - d)] . (3) 



The free energy of an elastic medium in a finite geom- 
etry can be computed by several different techniques for 
regularizing a divergent sum |p|,|6|,^|j2^] . The analysis re- 
duces to a ID problem for the in-planc Fourier modes 
^qi(^) = J dr±exp{-iq±r±)(l){r), where r_L = ix,y) 
and — [qx,qy). The partition function for each in- 
plane mode is given by 



du dv 'D(f)„ e 



[Wb.q^ + +5. 9i +9. v2] /fcfl T 



Here, Ti-b^q^ is the Fourier component of Tib and 
^ 'Dcj)^^ means that the functional integral should be 
taken over the paths that satisfy 



(4) 



(5) 
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and 



(<jO),<Jd) 



(6) 



where ' = d/dz. The path integral for a quadratic Hamil- 
tonian that contains a squared second-derivative has been 
calculated by Klcincrt [p6[ . Following his work, we de- 
compose each fluctuation path that satisfies the condi- 
tions (^) and (H) into two parts as <^q^ = 0m, + <50q^ , 
where (j)m.q^ = 0m. {z; u, v) is defined as the path that 
minimizes Ti-b.q^^ under the same boundary conditions. 
Then the Hamiltonian is decomposed as 



n. 



i/q I : 



d 

dz 





qj_ 



m,qj^ v^m, 



(9? 
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■■q^V^m.-q^ 



(7) 



(8) 
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The latter part can be expressed in terms of the bound- 
ary values as = a • M • a, where a = (ui, M2, f i, 1^2) 
and M = M{q^,d) is a 4 x 4 matrix whose calculation 
is straightforward. Accordingly, the partition function is 
factorized into "bulk" and "surface" parts as 



^q.= 

^b,qi = 



z 



s,qj_ ' 



,-Iq, //cbT 



0,0 



■^s,q^ = I dot 



-1/2 



Eu = diag(l, 1,0,0), 
Ev = diag(0, 0,1,1). 



(10) 
(11) 

(12) 

(13) 
(14) 



The bulk factor can be computed using Lagrange mul- 
tipliers that impose the boundary conditions Scjfq^ = 
Sd>' — IE™. In order to extract the interaction 



part of the free energy, we should regularize the par- 
tition hmction as Zq^{d) = liTaD^oci[Zqj^{d)Zq^{D — 
d)lZq^{Dl2)% so that ^q^(oo) = 1 (2|. Applying this 
to Zb^q^, we have 



(15) 




(16) 



In the same way, we obtain the regularized surface parti- 
tion function Zg^q^ , which has a voluminous expression 



that cannot be written down here. The interaction free 
energy per area is obtained as 



F : 



dqj 
{2n) 



In Zb. 



\nZ 



(17) 



First we focus on the bulk contribution by taking the 
"strong anchoring" limit ujs — + 00, — > 00, for which 
— 1. We show the result in terms of the disjoining 



^s,qj 



pressure 11 ~ —dF/dd. At the critical point po = 0, the 
contributions to the pressure from the regions q±^ < go 
and q±^ > go are given by, respectively. 



n< = - 



qod 



27rd3 
~ 2nd 
~ And 



ds 



s^(s — sin s cos s) 
— sm s 

qod < 1, 
qod > 1, 



(18) 



(19) 



where s is for \/%^-q^d, and 



n> = - 



27rd3 



dt 



t2(e-*sinhi + - t) 



sinh^ t — t"^ 



= -0.5463 ... X 



d3 ' 



(20) 



where t is for — q^d. Note that the total force 

n< + n> is proportional to \/d at large distances. To 
my knowledge, it is more long-ranged than any (ther- 
mal or quantum) Casimir interaction ever predicted for 
a 3D system. As seen from Eg. (|l8|) , the origin of the l/d 
tail lies in that every in-plane mode with a wavenumber 
qi^ < qo has a long-range correlation. It is further re- 
duced to the fact that the bulk soft modes of the system 
are distributed over the surface of the sphere |q| — go- 
This is in contrast to the case of a critical simple fluid, in 
which the soft mode is located on the single point q = 
and the Casimir force is proportional to d^^ (in the har- 
monic approximation). The extra factor d^ in the present 
system is thus explained by the dimensionality of distri- 
bution of soft modes in the q-space. Shown in Fig. 1(a) 
is the force profile at the critical point. 

In real systems, nonlinearity in the free energy makes 
the TL transition weakly first order pTf , and hence we 
do not expect to see the truly long-ranged force. For 
block copolymers, for instance, we estimate the typical 
and effective value of po/qo at the transition to be 0.4. 
Shown in Fig. 1(b) is the profile of the scaled pressure 
nd'^/fesT at finite values of po/qo- Even considerably 
above the transition, the force has a range of several 
times 7r/go, the lamella thickness. A weak oscillation 
arises from the region gj^ < go, while the contribution 
from the region gj^ > go is always monotonic. Note also 
the large magnitude of the interaction; the scaled pres- 
sure in its decay range is typically 10 times larger than 
the value Ild^ /ksT = —0.096 for nematic liquid crystals 
in the one-constant approximation 
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Next we study the role of fluctuations at boundaries. 
For the present model with a vanishing surface potential 
{u!s = 5s = 0), it is known that the mean-field solution 
^ = is destabilized above the bulk I-L transition 
The instability remains for a small but finite surface po- 
tential. I calculated the region of surface parameters 
for which the homogeneous mean-field solution is sta- 
ble at any temperature above the bulk I-L transition. 
It is plotted in Fig. 2(a) in terms of the dimensionless 
parameters fts = ujs/i^qo) ^^'^ ~ 9s/i^<lo)- Shown 
in Fig. 2(b) is the disjoining pressure for different values 
of (risjGs) in the stable region and at the bulk transi- 
tion temperature. The force has the asymptotic behavior 
lim^^oo n(i/(A;BTgg) = — C, where C is a positive num- 
ber that weakly depends on {ils,Gs)- An oscillation in 
n((i) appears for a small surface potential, and is ampli- 
fied as we approach the instability line. On the line, the 
quantity Yldj ikBTq^) is almost periodic in d except at 
small distances. The peak positions are approximately 
given hy d — {mr + 'ip)/qQ, where n is any integer and the 
phase shift ip moves over the region [0, tt] as we move over 
the instability line. To see the origin of the oscillation, we 
consider the simplifying limit (il^, Gg) = (oo, 0), which is 
an endpoint of the instability line. The interaction free 
energy in this limit reads 



F : 



2tt 



q±dq_i 



In 



2 sin 



+ In (l - e-V9i-9o'^) d{q^ - go) 



(21) 



Although it is a continuous function of d, its gradient 
(and so 11 (d)) logarithmically diverges at d = nTr/gg (> 
0), due to the contribution of the mode with {q±,qz) = 
{0,qo). This is closely related to the fact that the mode 
(j) (X sin((7o-z), which is soft in infinite bulk, is allowed 
by the boundary condition 0(0) = (/>((i) = only when 
d = riTr/qo- When we gradually increase d to cross one of 
the node points, the soft mode is abruptly allowed and 
then prohibited again, which corresponds to the rise and 
drop of the disjoining pressure. For physical (finite) val- 
ues of fls, boundary constraints are not strict and the 
oscillation in ll{d) is not singular [p8| . 

Up to now, we have assumed that the walls are ideally 
neutral. In general, the walls have a specific interaction 
with the fiuid, which is most simply described by the lin- 
ear Hamiltonian AHs — J dr hs(j) [S{z) + S{z — d)] . The 
surface field hg induces a mean-field interaction between 
the walls. Let us compare its magnitude with that of 
the Casimir interaction. To be specific, we consider sym- 
metric diblock copolymer melts with N — 10'^, a — 0.5 
nm, and T — 400 K, for which the lamella thickness 
Tr/qo equals 11 nm. For fig = Gg = hs/{eq^) = 1.0, 
Pq/qo = 0.4, and d — 4.77r/go = 50 nm, minimization 
of = Hb + Ti-s + ATis gives the surface value of the 
order parameter = 0.54 and the mean- field pressure 
n„if = —11 Pa |2^. On the other hand, the Casimir 
pressure for the same parameters and at the same dis- 



tance is Ilea = —69 Pa. The ratio Ilca/IImf = 6.3 shows 
that fluctuation effect on the total structural force is sig- 
niflcant within several lamella thicknesses. It should be 
mentioned that harmonic approximation to the free en- 
ergy is valid only for a weak surface fleld. It would be re- 
alized, for instance, by use of substrates coated by a thin 
sheet of random copolymers p^ . For a strong surface 
field, the order parameter is saturated at the walls where 
there appear wetting layers of the lamellar phase ||l^ . In 
that case, the above estimate should be applicable to the 
force between the surfaces of the wetting layers. 

To summarize, we predict an unconventionally strong 
thermal Casimir effect in structured fiuids above the 
isotropic-lamellar transition. At the ideally continuous 
transition, the force has a long-range tail proportional 
to l/d, which is due to the two-dimensional distribution 
of bulk soft modes in the q-space. The qualitative re- 
sult should commonly hold for model systems in which 
the bulk structure factor diverges at a finite wavenum- 
ber Although the force cannot be genuinely long- 

ranged in a real system, one should be able to access the 
temperature region where the force has a range of several 
layer thicknesses, within which it is much stronger than 
the conventional Casimir force at criticality. This implies 
that it also dominates other dispersion forces, such as the 
van der Waals interaction, which is often comparable to 
the critical Casimir force |0,||j9j. We have also studied 
the role of surface potential. Within a standard model 
of microemulsions, the force profile shows an oscillation 
near the surface-induced I-L transition. The oscillation 
originates in resonance of a soft mode in the slab geome- 
try. A measurement of structural force near the sponge- 
lamellar transition has reported an oscillatory force su- 
perimposed on an attractive background |3ll. However, 
more systematic study with controlled surface properties 
seems necessary to extract the fluctuation part of the in- 
teraction. I wish the present work to stimulate further 
experimental investigation. 
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FIG. 1. Profile of the Casimir force in the strong anchor- 
ing limit — > oo, gTs — > CX3. (a) Log- log plot of the absolute 
value |n| at the critical point po = 0. (b) Scaled disjoining 
pressure above the transition temperature. 



FIG. 2. (a) Region of stability of the mean-field solution 
(j) = Q. Outside the region, the surface-induced I-L transi- 
tion [15] preempts the bulk transition, even if there is only 
one wall (d oo). (b) Surface-potential dependence of 
the Casimir pressure. Plots are for {Q.s,Gs) ~ (0.5,0.5), 
(0.6,0.6), (1.5,1.5), and (cx),oo), all at the bulk transition 
temperature. The point (0.5, 0.5) is on the instability line. 
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